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FORTRAN PROGRAM FOR CALCULATING TOTAL-EFFICIENCY - SPECIFIC- 
SPEED CHARACTERISTICS OF CENTRIFUGAL COMPRESSORS 


by Michael R. Galvas 

Lewis Research Center and 
U. S. Army Air Mobility R&D Laboratory 

SUMMARY 

A method of predicting design point specific speed - efficiency characteristics of 
centrifugal compressors is presented with the computer program developed for the anal- 
ysis. The method is a one-dimensional mean streamline analysis conducted at fixed 
inlet stagnation conditions. Seven specific losses are calculated for each set of com- 
pressor geometric variables and inlet velocity diagram characteristics studied. The 
effect of these losses is then related to overall compressor performance and specific 
speed. By examining the program output the user can select values of inducer hub-tip 
diameter ratio, inducer tip-exit diameter ratio, impeller blade exit backsweep, impel- 
ler exit blade height-diameter ratio, and impeller exit absolute flow angle that will re- 
sult in maximum total efficiency for the chosen application. A numerical example is 
included to demonstrate compressor geometry selection for maximum efficiency. 


INTRODUCTION 

Compressor efficiency has been shown to be a function of specific speed. Specific 
speed is a characteristic that relates compressor inlet volume flow rate, rotation speed, 
and ideal enthalpy rise. High efficiencies are generally associated with high specific 
speeds and low efficiencies with low specific speeds. Compressor design point geom- 
etries that produce maximum attainable efficiency are also functions of specific speed. 
Knowing the variation in optimum design point geometric variables with specific speed 
permits rapid selection of high efficiency configurations. 

Balje (ref. 1) analytically correlated centrifugal compressor losses with specific 
speed and specific diameter. His major conclusions were that, for swirl-free centrifu- 
gal compressors, the optimum exit-inlet diameter ratio was primarily a function of spe- 
cific diameter and that backswept impellers produced highest efficiency in the specific 
speed range of 0. 70 to 1. 02. 



Another analytical method for predicting the variation in centrifugal compressor 
total efficiency with specific speed was described in reference 2. The velocity diagram 
characteristics and geometric variables that result in maximum total efficiency are pre- 
sented as functions of specific speed for several impeller tip speeds. A FORTRAN pro- 
gram was developed for the study reported in reference 2 and is given here with instruc- 
tions for its use. 

For given inlet stagnation conditions, the user can generate efficiency, pressure 
ratio, specific speed, and relative loss distribution data corresponding to various com- 
binations of impeller inlet velocity diagram characteristics and impeller overall geom- 
etries. By examining the output data a compressor geometry can be chosen which will 
yield maximum efficiency under the constraints imposed. The program can be used for 
working fluids other than air which approximate ideal gas behavior since the thermody- 
namic properties needed for the equations solved in the program are specified inputs. 


ANALYSIS 

The method of analysis is a one-dimensional mean streamline flow solution. Seven 
specific losses are calculated for each compressor configuration and specified inlet ve- 
locity diagram characteristics. These are inlet guide vane, blade loading, skin friction, 
disk friction, recirculation, vaneless diffuser, and vaned diffuser losses. Each of these 
individual losses is expressed as a decrement in compressor total efficiency. 

The enthalpy loss across the inlet guide vanes is computed at the rms inlet diameter 
using the equation for boundary layer losses presented in reference 3. 

The impeller losses due to blade loading and skin friction are calculated from the 
equations of reference 4. 

Impeller recirculation loss is computed using a modified form of the equation pre- 
sented in reference 4. 

Disk friction loss is calculated by the method of reference 5. 

Vaneless diffuser loss is determined by the numerical solution of the differential 
equations describing adiabatic flow in a radial passage derived in reference 6. These 
flow solutions were then used to solve the equation for total pressure presented in refer- 
ence 4. 

Vaned diffuser loss is calculated by determining pressure recovery attained in the 
vaned diffuser. Lines of maximum pressure recovery at a given area ratio were extrap- 
olated from test data reported in reference 7. The pressure recovery coefficient corre- 
sponding to an assumed exit Mach number of 0. 2 and throat conditions of Mach number 
and blockage was then determined by the iterative method described in reference 2. 


2 


INPUT INFORMATION 


Input information may be classified into the following different categories: (1) com- 
pressor geometry, (2) thermodynamic properties of the working fluid, (3) velocity dia- 
gram characteristics, and (4) iteration limits. The compressor geometry inputs are 

(1) inducer tip diameter, (2) inducer hub-tip diameter ratio, and (3) impeller exit back- 
sweep angles. The thermodynamic properties are (1) inlet stagnation temperature, 

(2) inlet stagnation pressure, (3) inlet stagnation dynamic viscosity, specific heat ratio, 
and gas constant of the working fluid, and (4) an estimated skin friction coefficient. The 
velocity diagram characteristics are (1) inducer tip speed, (2) inducer tip absolute crit- 
ical velocity ratio, (3) impeller exit-inlet tip relative velocity ratio, (4) inducer tip 
speed, and (5) prewhirl tangential-inducer velocity ratio. The prewhirl used in this 
analysis is solid-body vortex. Hereafter, "prewhirl velocity ratio" will be used in its 
discussion with the understanding that it is solid-body vortex. For iterations on inducer 
tip absolute critical velocity ratio the inducer tip speed is adjusted to preserve inlet ve- 
locity triangle similarity with that determined by the first pair of input inducer tip speed 
and inducer tip absolute critical velocity ratio. That is, the absolute and relative flow 
angles are held constant for successive iterations (see fig. 1). The input iteration limits 
are the numbers of values of (1) inducer tip absolute critical velocity ratios, (2) pre- 
whirl velocity ratios, (3) inducer tip-exit diameter ratios, (4) inducer hub-tip diameter 
ratios, and (5) impeller exit backsweep angles. 

A numerical example is included to demonstrate use of the FORTRAN program. For 
the example case, a compressor with a total pressure ratio of 6 to 1 and a mass flow 
rate of approximately 0. 9 kilogram per second was selected. Since high efficiency was 
of prime importance, a value of specific speed in the range of 0. 9 to 1. 0 was the prelim- 
inary target. This was determined by interpolation of the results reported in refer- 
ence 2. The target range of specific speed was used to select an inducer tip diameter of 
0. 0861 meter and rotative speed of 75 000 rpm. An inducer hub-tip diameter ratio of 



Figure 1. - Generalized inducer tip velocity triangle. 
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0. 3 was chosen on the basis of mechanical considerations. This represented a practical 
limit because of the large blade thicknesses encountered at the inducer hub. 

A one-dimensional continuity calculation (eqs. (B29) and (B31) to (B36)) indicated 
that an inlet Mach number of approximately 0. 7 (V 1T /V cr = 0. 628) would result in good 
inducer velocity diagrams with ample margin in mass flow to allow for blockage due to 
blade thickness. Inducer tip-exit diameter ratios in the range 0. 545 to 0. 565 were con- 
sidered for the pressure ratio of interest. Impeller exit backsweep angles of 0 to 45 
were studied to determine which backsweep would produce greatest efficiency at the 
6 to 1 pressure ratio. An inducer hub-tip diameter ratio of 0. 35 was studied to deter- 
mine whether aerodynamic improvements could be expected by increasing the selected 
ratio. 

Input information which corresponds to the sample case is given in the following 
table: 


Compressor geometry 

Inducer tip diameter, Djrp, 


0.0861 

.... 0. 3, 0. 35 


0. 54, 0. 545, 

0. 55, 0. 555, 0. 56 

Impeller exit backsweep angles, /? 2b , deg from meridional . . . 

... 0, 5, 10, 15, 20, 

25, 30, 35, 40, 45 

Thermodynamic properties 


Inlet stagnation pressure, Pq, N/m 


2 101 325.35 

Inlet stagnation temperature, Tq, 2 ®^ 15 

" 1. 788x10'° 

. . 1.4 


Inlet stagnation dynamic viscosity, pQ, (N)(sec)/m 
Specific heat ratio, y 


Skin friction coefficient, 0.004 

. 287.04 


Gas constant, R, J/(kg)(K) 


Velocity diagram characteristics 


Inducer tip absolute critical velocity ratio, V 1T' /V c 


Prewhirl velocity ratio, (V u /u)^ 

Impeller exit-inlet tip relative velocity ratio, W 2 /W 1T 
Inducer tip speed, m/sec 


0. 628 
. . . 0 
. . 0.7 
338. 14 


Iteration limits 


Number of values of inducer tip absolute critical velocity ratios 1 

Number of values of inducer tip-exit diameter ratios 5 

Number of values of inducer hub-tip diameter ratios 2 

Number of values of impeller exit backsweep angles 10 

Number of values of prewhirl velocity ratios * * 
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Figure 2. - Sample input sheet. 


Additionally, if a working fluid other than air is used in the analysis, an empirical 
equation expressing the dynamic viscosity of a function of temperature must be substi- 
tuted for those corresponding to equations (B59) and (B93). 

A sample input sheet is shown in figure 2. 


OUTPUT INFORMATION 


Before the computation of compressor performance is started, a complete list of 
input is printed out. Then, for each iteration the following information is tabulated: 

Compressor geometry 

Inducer tip- exit diameter ratio, D lT /D 2 
Inducer hub-tip diameter ratio, X 

Impeller exit backsweep angle, /3 2b , deg from meridional 

Impeller exit blade height-diameter ratio, b 2 /D 2 

Velocity diagram characteristics 

Inducer tip absolute critical velocity ratio V. /V 

J ? 1T / cr 

Prewhirl velocity ratio, (V u /u). 

Impeller exit absolute flow angle, £v 2 , deg from meridional 
Compressor performance characteristics 
Total efficiency, 77 
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Individual decrements in efficiency, A77 
Total pressure ratio, PR 
Specific speed, N g 
Head coefficient, if/ 

Sample output for the example case is shown in the next section. A survey of the 
output data for the example indicated that maximum efficiency at the approximate pres- 
sure ratio of interest was attained at an inducer tip-exit diameter ratio of 0. 555. Opti- 
mum impeller blade exit backsweep was 40° from meridional. Optimum impeller exit 
blade height -diameter ratio was 0.0306, and optimum inducer hub-tip diameter ratio was 
the initial selection of 0. 3. This combination of geometric variables resulted in a com- 
pressor total pressure ratio of 6. 07 and a total efficiency of 79. 1 percent. 

Head coefficient is an indication of the amount of deviation from maximum total ef- 
ficiency at a given specific speed. By comparing the calculated value to a value inter- 
polated from the figures of optimum head coefficient as a function of specific speed from 
reference 2, the relative penalty in off -optimum efficiency can be estimated. 


SAMPLE OUTPUT 


ofcfciN fcXECOT ICR* 
GAM 

1*4000 
CF 
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ERROR MESSAGES 


Two major errors can be encountered in the use of this program. A poor combina 
tion of inlet flow velocity, prewhirl velocity ratio, inducer tip speed, impeller inlet tip 
exit diameter ratio, and impeller relative velocity ratio will result in impeller exit 
velocity diagram characteristics for which a velocity triangle cannot be calculated. In 
this case calculation of the flow solution downstream of the impeller exit is deleted and 
the printout for the irrational combination of variables is suppressed. The second major 
error is the lack of convergence of the calculated values of vaned diffuser area ratio re 
quired to decelerate the flow to an exit Mach number of 0. 2. Values of vaned diffuser 
inlet blockage and Mach number can result in the interpolation of pressure recovery 
coefficients for which the specified exit Mach number cannot be attained within an area 
ratio of 5. 0. When this happens, AREA RATIO IS NOT IN BOUNDS is printed out. 


FORTRAN PROGRAM 

The FORTRAN listing, input and output samples, and definitions of the FORTRAN 
variables are given in this section. For each combination of input variables the main 
program calculates overall compressor design point performance from empirical loss 
estimates and adiabatic flow relations using a mean streamline one-dimensional anal- 
ysis. 


Program Input 

Input variables are shown in this section with the definitions of the FORTRAN van- 
ables used. These variables appear in the section Main Program FORTRAN Variables 
and Engineering Symbols but are repeated here for the convenience of the user. 

GAM specific heat ratio, y 

2 

POP inlet stagnation pressure, Pq, N/m 

POP inlet stagnation temperature, Tq, K 

rGAS gas constant, R, J/(kg)(K) 

UlT inducer tip speed, Upp, m/sec 

D1T inducer tip diameter, D-^, m 

jyjUO inlet stagnation dynamic viscosity, \i q 

CP skin friction coefficient, Cj 
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VOVCR(I) 

VUT(J) 

DRAT(K) 

LAM(L) 

B2(M) 

NVOVCR 

NVUT 

NLAM 

NB2 

W20W1T 


inducer tip absolute critical velocity ratio, V 1T / V cr 

inlet solid body prewhirl tangential- inducer velocity ratio, (V u /u)j 

inlet tip-exit diameter ratio, D^rp/Dg 

hub-tip diameter ratio, X 

impeller exit blade angle, 

number of values of Vjrp/V cr of interest 

number of values of (V u /u) ^ of interest 

number of values of X of interest 

number of values of of interest 

impeller exit-inlet tip relative velocity ratio, Wg/Wjrp 


FORTRAN Listing 


*IBF TC MAIN 

C FORTRAN PRCC-RAM FOR CALCULATING TOTAL EFFICIENCY-SPECIFIC SPEEO 

C CHARACTERISTICS OF CENTRIFUGAL COMPRESSORS 
DIMENSION VC VCR (Qi , VL 1 1 8 ) 

DIMENSION DR AT I 16) ,LAM ( 16 ) , BE TA2 ( 16 ) 

DIMENSION P3F(20) , RARRAY ( 2C ) , XMARR ( 20 ) 

CIMENS ION FI 10) ,S<20) 

DIMENSION AMT(4),BARR(6) , PR EC 1(4, 6) ,PREC2(4,6) ,PREC3(4,6) 
DIMENSION PREC414»6) » PREC5 (4 » 6 ) 

EXTERNAL FCNCT 
REAL LAM, MUO »LAMX»LCG 
1 R EAC 1 5 ,510 ) OAM»RGAS»PCF,TCP»U1T»D1T,MUO 
REACI5 ,5 10 ) CF > N20nlT 

REACI5 ,511 ) N VU T, NLAM , NCR AT , NVOVCR ,N82 
REACI5 »51C ) (DRAT( I ), I=1,NDRAT) 

REAC (5,510) (LAMII), 1*1, NLAM) 

REACI5 »51G ) (BETA2I I ) ,i=l,NB2) 

R EACI 5 ,5 1C ) IVUTII), 1*1, NVUT) 

READ! 5,510) ( V0 VCR ( I ) ,1*1, NVOVCR) 

NR IT El 6, 520) GAM,KGAS,POP,TOP,C1T,D1T,MUC 

WRITE I 6, 521 ) CF, W20w 1T,NVUT,NL AM, NC RAT, NVOVCR, NB2 

NR I TE ( 6, 522 ) (DRAT (I) , 1*1 ,N0R AT ) 

mR ITE ( 6 , 522 ) ( LAM( I ) , 1=1 , NL AM ) 

wRITE(6,524) (BETA2I I),I=1,N62) 

NRITE(6,525) I VLT( I ) , 1=1, NVUT ) 

WRITE! 6,526) <VCVCR( I ) , 1*1, NVOVCR) 

DATA! AMT ( 1 ) , I* 1 ,4) / . 2 , .4 , . 6 , . 8 / 

DATA(8ARR( I) , I* 1 , 6 ) / ,02, .04 , .06 , .0 8 ,» 10 , . 12/ 

DAT A 1 1 PREC 1 1 1, J ),I=1, 4) , J *1 , 6 ) /. 234 , . 244 , .257 , .269 ,. 215 , . 2 24 , . 233 
1 • 24 3 1 • 20 7 t . 2 25 f • 22 3 » • 232 m193m 199 r*206f*212 ? » 1 83 t *190#*196f* 2G 2 » 
2169 9 • ! 76 » • 162 9 • 168/ 

CATA( ( PREC2 ( 1 1 J > t 1 = 1 9 4 ) f J = 1 , 6 i /.644, • 6 70 9 *696 9 .722 » .620 9 .638 1 • 656 



1.674.. 590..6C6,.623» .£39 « . 562 » .576 ,. 590 , . 6C5 ,.538 , . 551 » . 564 , . 578 >. 

2510.. 524..53fc,.55 2/ 

DATA! ( PREC3 l l* J ), I =1 , 4 ) , J=1 ,6 ) /. 782, • 785, .79 6 , .802 , . 750 , . 75 6, • 762, 
1.768*. 70 8*. 716, .724 » » 132 1. 672* .680 *.687 *.695 *.652 * . 648 * . 654 * • 660 * . 
260 4*. 612 *.615**626/ 

U AT A ( ( PKEC4 I I * J I* 1 = 1*4) *3=1*61/. 842 * . 838 * . 83 3*. 828*. 8*. 8*. 8*. 8*. 75 
12*. 75 6 * . 760* .763* .7 10 *.713*. 716*. 719*. 675*. 678*. 680*. 683*. 6 30 *.635 
2**6 40 * .646/ 

C A T A ( I PR EC 5 1 1»J I* 1=1*4) *J = 1*6) /.676». 66 5**852* .838 *.832**825* .618* 
1.612*. 78, .78 ,.7 £,.78*. 736*. 735*. 735*. 754*. 652*. 654*. 695*. 696*. 644, 
2. 647*. 6 5b* .652/ 

02 = GAM-1 • 

Ol=GAM*l. 

CF=GAM*RGAS/C2 
DC 100 I = 1 * N VG VCR 
DC 103 J = 1 * N VuT 

CMEGA=2.*U1T*V0VCR( I ) /VC VCR ( 1 ) /D1T 
VCR=$3RT(2.*GAM/G1*RGAS*TCP) 

V 1 1 =VC VCR ( 1 } *VCR 

AL1=ARS1N(VLT( J)*01T*VGVCRI I ) / VC VCR (1 ) / VlT ) 

S IN A=S INIAL1 ) 

COS A=CCS I AL1 ) 

VM1 =V1 T*CG 5A 

RC17 = U1T*11.-VCT( 3) )*VGVCR( I ) / VC VCR ID 

ol=ATAN(fcblT/VMl) 

wlT = VMl**2 + VHT**2 

n1T = S3RTU1T ) 

AKJ=VM1**2*2 .*< VUT I 3)*LlT*VCVCRt I ) /VOVCRI 1) ) **2 
w2* W20R1T IT 
CC ICO K=1 * N CRA T 
WRITE 1 6*500 ) 

WRITE 16, 501) V0TI3) , VCVCR (I) , DRAT ( K ) 

0RT = 1./CRATU<) 

u2=ClT*VCVCR(l)/V0VCRCi)*CRT 

02 = C1TA0RT 

OC ICO L = 1 *N LAM 

wR I TE I 6* 502) 

wR ITEI 6* 503 ) LAP ( L ) 

wR I T£ 1 6 * 5C4 ) 

LAMX=LAM(L ) 

01H=LAMX*C1T<VCVCk{ 1 ) /VCVCR I 1 ) 

D1Fi=LAMX*C1T 

VC1H = LAMX*VCT( J)*U1T* VCVCR ( I )/ VO VCR III 
VM1H=3CRTI XKJ-2.*VClh**2) 

Rllh=blH-VUiH 

blh=ATANUClF/VMlH) 

Wlh-S3RTI VMlF**2+wUlF**2) 

DMF=SURT{ClT**2*il.*LAMX**2 )/2.) 

UlMF=01T*VGVCR( II/VCVCR(1)*DMF/D1T 
VU1MF=0MF/C1T*VUT( J ) *C1T*VCVCR i I )/ VCVCR (1) 
VMlMF=SgRTl*Kj-2.*VClPF**2) 

AL1MF*ATAM VClMF/VMlKF) 

V1MF=SwRTIVC1MF**2+VM1MF**2) 

WC1MF=C1MF-VL1MF 
wlMF=SCRTl VMlMF**2+hClMF**2) 
bIMF=ATANUClPF/VMiMF ) 

TlM=T0P-VlMF**2/2./CF 
E1AV= 1 61 *d 1 PF *8 1H ) / 2 . 

Dh I GV=0. 
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P1P=PCP 

PPCPl=(l.*VlPF**2/2./CP/TlM)**(GAM/G2) 

PCPP1=1./PPCF1 

P1=P1P*PCPP1 

R1P=P1P/RGAS/T0P 

R1 = R1P*(P1/P1PJ**(1./GAM) 

0=2. 1* 159*011**2* U.-LAMX**2H‘VMlMF/4. 
Sh=o*Rl 

RO P=POP/RG AS / TCP 
RE=U2*C2/MCC*RCP 
TlPP=TlM*hlMF**2/2./CP 


T2PP=TiPP+(C2**2-UlMF**2 )/2./CP 
T2=T2PP— R2**2/2./CP 


A2=SORTl GAM*FGAS*T2) 

PHl=VrtiMF/C2 

R0C2=( PHI**2*(DRF/D2 ) **2*R2CR 1 T**2* ( PH l**2+DRAT ( K > **2) ) /2. 

IF(VUT(J}.LT..01) GC TC SO 

ALSTAG=ALlMF/2. 

SINX=S1N(AL1PF) 

CGSX=CC$IAL1PF) 

fcS=.CO 76/(CC5X-.C25]*(1.*CCS(ALSTAG)/.7) 

AKE=V1MF**2/ 2. 

AKEIL=AKE/(1 «— E S * 

PGPP1=(1.-AKE/CF/TCP)**( GAM/u2 J 
P1GPP0=( l.-AKEID/CP/TOPJ **(GAM/G2J 
P1P=P0P*F1CFFC/PCPP1 
R1P=P1P/RGAS/TGP 


P1=P1P*PCPP1 

HI = R1P*(P1/P1P)**(1./GAM) 

sw=o«ri 

RV=R1*VM1MF 

PArA=RWROP/VCR 

AMSTAR=G.1G 

li Y=FUNCT( AflSTAR,Gl,G2)-PARA 
AMSTAR=AMSTAR*.GOi 
1F4Y.LT.C.) GO TC 11 
VC=AMSTA«*VCR 


AMC=9.79t5E-7*TlM**1.5/( T1P+11C.A) 

R EC = S W /AMU /C IT 

DHIGV=C.4*S I AX* ( VO**2 ♦ V1KF**2 i /2./COSX/REO**.2 
9C CONTINUE 

DC 1CU M=1,NE2 
62X=6ETA2(M)*.C1745 


Z=fe.5*( l.*CMF/02i/( 1.-CRF/C21* 
BPSLIM=l./EXP(8.i6*CCS(B2X) /Z i 
VSL = SgRT (CCS (S2XJ )*02/Z**.7 
i F ( (CMF/C2) • GT . EP SL I H ) VSL=U2* 
lEPSLXMj/J l.-EPSLIH) )**2 ) *U2* ( ( 
IF((VSL*CGS(E2X)/w2i.G7.1.) GO 


CGS ( (61AV*62X)/2.» 


(SgRT (COS (62X1 )/Z**.7)*(l. 

CMF/D2— EPSLI P) /(1.-EPSLIM) 
TC 1C1 


-( < CMF/D2- 
) **3 


D£LTA=ARS I A ( VSL *CCS 1 B2X ) /*2 i 
rt2ID=0 2*COS(fi2X*DELTA)/CQS< B2X) 
Vu2=U2-VSL-R2ID*SIMb2X) 


VM2=r2IL*CCS(B2X) 

IF( VR2.LT.C. J GG TO 1C1 
IF(VU2.LT.O. 1 GC TC 1C1 
AL2=AT AN( \rU2/VM2i 
RU2=U2— VC2 


V 2 = SORT ( WU2**2*VM2**2 ) 
T2P=T2*V2**2/2./CP 



GhAERG=CP*TGF*l T2P/TCP-1. ) 
t - TH=LHAERG/L2**2 

UF* «aT/G2*<Z/3.1415S*(l.-GlT/02)+2.*GlT/D2 > 
GF=1.-*2C*1T ♦C.75*gTH/OF 
0H6L=0 .C5*CF**2*G2**2 
0HRC=0 .02*StikT ( TAN ( AL2 ) ) *DF**2*U2**2 

LuO*( l.-CMF/ •3048)/CCStB2X)/2» 

R2G=RA*< T2/T 1M)**1 1./G2) 

5C RhG2*R2G 

BBC F* «C13 56 l*RHQ2*u2**3*Q2**2 /Sk/RE**»2 


B2 = Sk</ ( 3« 141 59*RH02*C 2*VR2 J 
GHYD=Z/3.14l59/CCS<B2*)+D2/e2 
DHYC*l./CHYG*DlT/02/(2./ll.-LAMX) 
1*SGRT( l.*(l.*LAMX**2) /2.M1ANIB1) 


♦2,*Z/3.14159/(1.*LAMXI 
)**2) I 


0HSF=5.6*CF*LGD/DHYD*hGG2*L2**2 
GhACT = CHAERG«OHCF«-CHRC 
GHiO*UhA£RO-LHI GV-GhBL-GHSF 


5i 


ETAR=UhIC/ChAERC 

TX=ETAK*QhAEFC/CP/TOP*l. 

P2P=Tx**( GAM/G2 J*P1P 


P2=P2P*( T2P/T2)**1-GAK/G2) 
R2G*P2/RGAi / 12 

IF( A6SUR2G-FHU2)/R2G).GT«.CC1) GG 


TO 


5C 


XM2 = V2/A2 
R2=C2/c. 

AML*S.7965E-7*T2**1.5/IT2+11C.4J 

ANL*AML/R2G 

BO=l. 

XM=XM2 
ALPhA=AL2 
NC = 1 
k*l.Q 

F ( i I »XM2** 3/ (l.*G2/2.*XM2**2J **(GAM/G2I*R 

P3P 1 1 ) -P2P 
RARRAY (1 ) =1 «C 


XMARRi 1)=XM2 
GELTAR=C»02 
ZETA*CF*H2/B2 
2C NC=NG*1 
JAN* i 


XM1*XM 
ALPH1 =ALPhA 

21 QELTAS=C«37*CELTAR*R2/8. /CG S C ALPHA i /( V2*DELT AR*R2/ CCS ( ALPHA )/ANU>* 
1* • 2 

6=eC-2.*CELTAS/B2 

GtLTAB=bO-B 

VARM*-2«*( 1« ♦G2/2»*XM**2 )/(XM**2-l»/CGS(ALPHA)**2) *( ( GAM*XM**2-T AN 
It ALPHA 1**2 i *ZETA/BO /CCS (ALPHA )-!• /BC*GELT AB /CELT AR-1 ./COS ( ALPHA ) ** 
22/R)*XM**2*CELTAR 

VARAL=1./CCS(ALPHA)**2/1XM**2-(1./C0S(ALPHA» )**2I* ( < 1. +G2/2. *XM**2 
l)*ZETA/bG/CGS( ALPHA )-1./B0 *CELTAd/GElTAR-XH**2/R>*TAN( ALPHA »*DELTA 

2R 

BC = B 

IF(JAN.bw.l) VARAL1=VARAL 
IFl JAN. EL. 1) VARMi=VARM 
lF( jAN.tG.i ) GC TO 22 
xh=xm**2+varf 

XP = S%RT t XH ) 

TAlPH=TAN( ALPHA )»VARAL 
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AlPHA*ATAN(TALPH) 

JAN=JAN+1 
K=R*DtLTAR 
GO TC 21 

22 VARAL=(VARALl*VARAL)/2. 
TALPF=TAN(ALFH1 )+VARA L 
ALPHA=ATAN(TALPH) 

VARM=( VARMl*VARM)/2. 

XM=XMi**c*VARM 

XM»SW*T(XM) 

ACCCSK=l./( l.+G2/2.*XM**2 ) 


RHCR=l./( l.*G2/2.*XM**2 )**(1./G2 ) 
F (NC) =XM**3*ACQLSR*RFCR*R 


IF ( NO .fcg.2 ) S( NO) = ( H NO) *F ( NC-1) i *C. 5* CEL TAR 
IF(N0.Eg.2) GO TO 23 
CALL FNTGRKNu, .02. F, S) 

23 TPL=1./(1 .♦GAM*CF*R2*S(N0 ) /CC S ( AL2 ) /B2/ XM2* ( 1. ♦G2/2,*XM2**2 )**(6AM 
1/G2) ) 


PP*TPL*P2P 
P3P (NU)=PP 


RARRAY (NC )=R 
XMARR(NG )=XM 
IF( XM.GT.C.fi J GO TC 2C 

RVANE=< XMARR (NO-1 J -0 . 8 ) / ( XMARR ( NL-1 ) -XMARR (NC) I *( RARRAY ( NO ) -RARRAY 
1( NC— 1) j *RARR AY ( NO-1 J 

PTHP= ( R VAN E- RARRAY ( NC- 1 ) )/.02*(P3P( NO— 1 )-P3P(NC) ) 
PThP=P3P(NG-l)-PTHP 

XMACH=(RVANE-R ARRAY (NC-1) ) /.C 2*( XMARR ( NO-1 )- XMARR ( NC ) ) 

XMACH-XMARR ( NO- 1 l-XMACH 

IF ( XMARR ( 1 ) .LE.C.8) XMACF=XMAKR ( 2 ) 

IF(XMARR(1).LE.0.8) PTHP=P3P(b) 

PTH=PTHP/( 1. ♦G2/2.*XMACh**2)**(GAM/G2) 

uHVLC=CP*T2P*( ( PTH/PThP)**(G2/GAM)-(PTh/P2P)**(G2/GAM) ) 

eT=l.-B 

AR=XHACH/C.2*( ( l.«-G2/2.*C.C4)/(l.*G2/2.*XMACF**2))**(Gi/2./G2) 
PP£XIT=PThP 
7 l AR1=AR 

ARNUM*AR*PPEXIT 

IF( (AK— 1.2)«CT*0..ANC#(AR-2.).LT.fc.) GO TO 60 
IF( (AR-2. ) .GT.O..ANG. (AR-3. I.LT.C. ) GO TO 61 
IF( (AR-3.). GT.O.. AND. (AR-4 .).lT.G.) GO TC 62 
IF( ( AR-4.).G1.C..ANG. (AR-5. ).LT.G. ) GO TC 63 
mRIT£(6,11C) 

11C FORMAT ( 28h A FEA RATIO IS NOT IN BCCNDS) 

GO TO 102 

6C CALL LININT (XM ACH , BT , AF T ,6 A PR , PREC1 , 4 , 6 , FI ) 

CALL LININT (XMACh , BT ,AMT ,BARR .PREC2 , 4, 6, F2) 

CPSTAR = ( AR-1 ,2)/.8*(F2-F 1 ) *F 1 
GO TO 6B 

61 CALL LININT ( XM ACh « BT . AMT .BARR .PREC2 . 4 . 6 ,F 2 ) 

CALL lININT (XMACH#BT,AMT,£ARR|PR£C3»4,6,F3) 

CPSTAR=( AR-2.)* (F3-F2) +F2 

GC TO 68 

62 CALL LININT (XM ACH » BT ,AMT , BARR , PREC3 » 4 , 6 ,F3 ) 

CALL LININT (XMACH, BT, AMT , d ARR ,PRtC4 , 4 , 6 , F4) 

CPSTAR=(AR-2.)*(F4— F3)+F3 

GC TC 68 

63 CALL LININT (XMACh , BT ,AMT , BARR , PREC4 , 4 , 6 ,F4 ) 

CALL LININT ( XM ACh $ BT , AM T » B A RR , PRt C5 , 4. 6 , F 5) 

CPSTAR*(AR-4.)*(F5-F4 )*F4 


6b PEXiT=CPSTAK*l PThP-PTh)+PTH 

PPEX IT-PEX i T * t l.+G2/2.*t.C4 )** t GAM/G2) 

AR=ARNCM/PPEXIT 

iF(A6MARl-AR).LT..CCl) GO 70 69 

69 DHVO C =T2P*CP*( (PEXIT/PPEXIT)**IG2/GAM)-(PEXIT/PTFIP )**(G2/GAMJ ) 
cTAO=( GHAERC-DFlSF-CHBL-DhVLD-DFiVD-UHlGV) /DHACT 


B2C2=B2/C2 
PR=PP£XIT /POP 
FliU = ETAC*£HAERO 
PM=Hl C/L2**2 

XN5=CMEGA*SCRT (C)/HIC**.75 

U£ ItV-CH l GV/CHACT 

UttJL = OFBL/CFACT 

DtSF=UhSF/ChACT 

OECF=CFDF/CFACT 

DERC-GFRC/CFACT 

DEVLC = tFiVLC / CHACT 

DEVu=UhVC/CF ACT 

WRITE! 6?5C5)^BETA21M) ,C£ IGV ,DE BL , DESF ,D£OF ,DERC»CE VLD , DEVD , PR ,XNS» 

5i,t 1 |0RM;T( 2 ,P !3fSlNLET PREkFIRL C VO/OU .5X.51HI PPELLER INLET TIP ABSO 
iLOTE CRITICAL VELOCITY RAT IC , 5X» 34HI MPELLER INLET-EXIT DIAMETER RA 
2TILI 

501 FCRMAT(1CX,F1C.4,32X,F1C.4,32X,F10.4) 

502 FORMAT (39X,29Fi INLET hUB-T IP DIAMETER RATIO) 

5 C -£ FORMAT (48X » F 10. 4) ccc ir itwr v rFrfiFMFNTS* 

50A FORMAT 19H IMF ELLER , 2X ,54H EFFICIENCY CECREMENT5 

, 3 x . BhPRESSURE » 3X » BHBPEC I F IC »25X» AHREAD »9X» 5HT OTAL/ » l 

20 H cACKSrEEP »3X »3HIG V , 5X* 2F.BL , 6X ,2HSF , 6X , 2HDF , 6X ,2HRC ,6X,3FIVLC * 6X» 
32HVD»6X»5FiRA1IC»7X»5FSPEED»4Xf 5HB2/D2 , 4X.7HALPHA 2 , 2X , 11HCGEFF IC I 
4EM,2X,1QFEFFICIENCY) , , v £i t 5< c . , 

5C5 FCRMATIF1CS.4,1X,F6.4,2X,F6.4»2X,F6.4,2X,F6.4»2X,F6.4,2X,F6.4,2X,F6 

1.4, 3X,F3.3,3X,F 8. 4 , 3X ,F 6 . 4 ,3X , F7. 3 ,2X ,F10.4, 5X,F6.4) 

5 1C FORMAT ( 8 F 10 • M 


520 FORMAT 14X^3FGAM t 12X ,4hRGAS »12X,3HPCP ,12X, 3HTCP ,12X » 3HU1T , 12X , 3HD1T 

1, 12X 1 3FIMUC /*1X,F1G.4,5X»F1C.4,5X,F1C.3,5X,F1C.4,5X,F10.4,5X,F1O.4» 

521^F0RMAT I4X, 2FCF , 12X ,6Hk2Cii»lT , 10X,4HNVUT ,llX»4hNLAM , 11X, 5HN0RAT , 10 X, 
16HN VGVCR , 1CX ,3HNB2/ , 1X,F 10.4 , 5X ,F1C. 4, 8X , 14 , 11X , 14 , 12X , 14,1 IX, 14,1 

DRAT ARRAY / 5X,8F10.4/5X,8F10.4 ) 

LAM ARRAY/5X,eFi0.4/5X,BFlC.4) 

B£TA2 ARRAY/5X, 8F 10.4/ , 5X , 8F 10 .4 ) 

WLT ARRAY/5X,8F1C,4) 

VGVCR ARRAY /5X,8F1C.4I 


5^.2 

52j 

524 

525 

526 
1C2 
101 
lot 


FORMAT (IX, 11 F 
FORMAT ( IX , 10 F 
FORMAT I IX , 12 F 
FORMAT liX, 1CF 
FORMAT 1 1X,12F 
COM INOE 
COM I NOE 
CENT INOE 


Gii TC X 
ENC 
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Function FUNCKAMSTAR, Gl, G2) 


This routine calculates (pV/p'V cr )Q for an estimated Vq upstream of the impeller 
inlet guide vanes. The solution obtained by trial-and -error process is used to compute 
the total pressure loss across the inlet guide vanes. 

FUNCT (pV/p’V cr ) 0 

AMSTAR V Q /V cr 
Gl y + 1 

G2 y - 1 

JIbFTC FUNCT 

FUNCTION FLKCTI AMSTAR,tl,G2> 

FUNCT* ( l»-G2/Gl*AMSTAR**2)**( 1 » /G2 ) * AMS TAR 

RETURN 

tNL 


S ubroutine UNINTCXl, Yl, X, Y, TN, MX, MY, F) 

This subroutine interpolates a value of maximum pressure recovery coefficient C** 
from a table of vaned diffuser throat Mach numbers and aerodynamic blockage given as 
input. 

XI input, vaned diffuser throat Mach number 
Yl input, throat aerodynamic blockage 

X input array of throat Mach numbers 

Y input array of throat blockages 

TN input two-dimensional array of C** corresponding to throat Mach numbers and 
blockages 

MX input, number of throat Mach numbers 
MY input, number of throat blockages 
F output, interpolated value of C** 
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UaFIC LININT 


GOtiKCcTlNE L lNIMtXit Y1,X,Y ,TN»MX, ►Y,FJ 

DIMENSION X(PA) ,YiMY) »TMMX,MYJ 
DC 10 J3*2»PX 

1C IFUl.LE.Xl J2)i GO TO EC 
J3=KX 

20 DC 30 J4*2.PY 
30 IFtYi.Lt.Yi JO) GO TC 40 
J4*KY 

4C J1=J3— 1 
j2=J4-1 

EPSl=(Xl-X( Ji>)/(X1J2)-X( Jl») 

£PS2=l Yl-Yl Ji) ) /( Y( J4 )-Y< J2 ) ) 

tPSi=i.-EPSl 

EPS4=i.-EPS2 

F=TN(JJ*J2I*EPS3*EPS4+TN(J3»J2)*EPS1*EPS4+TN<J1»J4 ) *EPS2*EPS3* 
1TM J3, J4)*fcPSl*EPi>2 
RETORN 
ENC 


Subroutine FNTGRUNO,. 02, F, S) 

This integration subprogram from the IBM scientific subroutine package is used to 
integrate the function M 3 (pa/p'a') in the vaneless diffuser. 

NO number of equally spaced radii 
. 02 radius ratio between stations 
F M 3 (pa/p'a’) 

S integral of M 3 (pa/p'a') 


Main Program FORTRAN Variables and Engineering Symbols 


ACOUSR 

a/a' 

AKE 

KE 

AKEDD 

^id 

ALPHA 

a inside vane less 

ALPH1 

temporary storage 

ALSTAG 

a st 

AL1 

ff lT 

AL1MF 

^lMF 


diffuser 
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AL2 

a 2 

AMSTAR 

M* 

AMT 

array of throat Mach numbers 

AMU 

M 

ANU 

V 

AR 

AR 

ARNUM 

ARp^ 

AR1 

temporary storage 

A2 

a 2 

B 

B 

BT 

vaned diffuser inlet blockage 

BARR 

array of vaned diffuser inlet blockage 

BETA2 

/3 2b array 

BO 

B 0 

B1 

% T 

B1AV 

^1AV 

B1H 

^1H 

B1MF 

^1MF 

B2 

b 2 

B2D2 

V D 2 

B2X 

^2b 

CF 

C f 

COSA 

COS 

COSX 

cos a lMF 

CP 

C 


P 

CPSTAR 

c r 

DEBL 

Ar ?BL 

DEDF 

a? ?df 

DEIGV 

a ^igv 

DELTA 

6 


DELTAB 

AB 

deltar 

AR 

DELTAS 

A 6* 

DERC 

a?? rc 

DESF 

At? sf 

DEVD 

At ?vd 

DEVLD 

At ?vld 

DF 

D f 

DHACT 

Ah act 

DHAERO 

Ah aero 

DHBL 

Ah BL 

DHDF 

Ah DF 

DH1D 

AH id 

DHIGV 

Ah IGV 

DHRC 

Ah RC 

DHSF 

Ah SF 

DHVD 

Ah VD 

DHVLD 

Ah VLD 

DHYD 

D hyd/ D 2 

DMF 

°MF 

DRAT 

^1T^2 

DRT 

D 2 /Dit 

D1H 

d ih 

D 1T 

d it 

D2 

D 2 

EPSLIM 

€ lim 

ES 

e s 

ETAD 

T ?AD 

ETAR 

^R 

F 

M 3 (a/a' 


22 


FI 

C** at vaned diffuser area ratio of 1. 2 

F2 

C** at vaned diffuser area ratio of 2. 0 

F3 

C** at vaned diffuser area ratio of 3. 0 

F4 

C** at vaned diffuser area ratio of 4. 0 

F5 

C** at vaned diffuser area ratio of 5. 0 

GAM 

r 

G1 

y + 1 

G2 

y - 1 

HID 

H id 

JAN 

iteration counter 

LAM 

X array 

LAMX 

X 

LOD 

L/D 2 

MUO 


NB2 

iteration limit 

NDRAT 

iteration limit 

NLAM 

iteration limit 

NO 

increment counter 

NVOVCR 

iteration limit 

NVUT 

iteration limit 

OMEGA 

O) 

PARA 

(pv/p’V cr ) 0 

PEXIT 

p 4 

PHI 

<p 

POPP1 

(P/P’)! 

PP 

P’ 

PPEXIT 

?4 

PPOP1 

(P’/P)! 

PR 

P4/P0 

PREC1 

array of C** at area ratio 1. 2 


23 



PREC2 

array of C** 

at area ratio 2. 0 

PREC3 

array of C** 

at area ratio 3. 0 

PREC4 

array of C** 

at area ratio 4. 0 

PREC5 

array of C» 

at area ratio 5. 0 

PSI 

* 


PTH 

P 3 


PTHP 

P 3 


POP 

Po 


PI 

P 1MF 


PIOPPO 

Pimf/po 


PIP 

P 1MF 


P2 

P2 


P2P 

„ T 
p 2 


P3P 

_ T 

p 3 


Q 

Q 


QTH 

^th 


R 

R 



RARRAY temporary storage 
RE Re 

REO Re based on inlet stagnation conditions 

RGAS R 

RHOR p/p' 

RH02 p 2 

RV ipV) 

RVANE r 3 

ROP p' 0 

Ri P 1MF 

R1P P 1MF 

R2 r 2 

R2G p 2est 
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s 

integral of M 3 (pa/p'a T ) 

SINA 

sin a 1T 

SINX 

Sin “IMF 

SW 

w 

TALPH 

tan a 

TPL 

P’/Pfe 

TX 

^R Ah aero/ C p T 0 + 1 

TOP 

T b 

TIM 

T 

IMF 

T1PP 

rp » T 
A 1 

T2 

T 

2 

T2P 

T' 

2 

T2PP 

rn T * 
2 

U1H 

U 1H 

U1MF 

U 1MF 

U1T 

U 1T 

U2 

u 2 

VARAL 

A tan a 

VARA LI 

temporary storage 

VARM 

AM 2 

VARM1 

temporary storage 

VCR 

V cr 

VM1 

V mlT 

VM1H 

V mlH 

VM1MF 

V mlMF 

VM2 

V m2 

VOVCR 

V lT/ y cr ^ay 

VSL 

V SL 

VUT 

(V u /u)i array 

VU1H 

V ulH 

VU1MF 

V ulMF 

VU2 

CM 

VO 

v o 
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V1MF 

V 1MF 

V1T 

V 1T 

V2 

V 2 

WOU2 

< w /“2>av 

WTJlH 

W ulH 

WU1MF 

W ulMF 

WU1T 

W ulT 

WU2 

W u2 

W1H 

W 1H 

WlMF 

W 1MF 

WIT 

Wit 

W2 

W 2 

W2ID 

W 2id 

W20W1T 

w 2 /w 1T 

XKJ 

K 

XM 

M 

XM1 

temporary storage 

XM2 

m 2 

XMACH 

m 3 

XMARR 

array of vaneless c 

XNS 

N s 

Y 

dummy variable 

Z 

Z 

ZETA 

C 


M's 


Lewis Research Center, 

National Aeronautics and Space Administration, 
and 

U. S. Army Air Mobility R&D Laboratory, 
Cleveland, Ohio, April 18, 1972, 

132-15. 
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APPENDIX A 


AR 


a 

B 



H 

Ah 

K 


KE 

L 

M 

N s 

PR 


P 

Q 

q 

R 

R 

Re 


r 

T 


SYMBOLS 


area ratio 

local acoustic velocity, m/sec 

diffuser effective depth ratio 

vaned diffuser throat aerodynamic blockage 

blade height, m 

skin friction coefficient 

specific heat at constant pressure, J/(kg)(K) 

maximum pressure recovery coefficient at a given area ratio 
diameter, m 
diffusion factor 

inlet guide vane loss coefficient 

overall compressor enthalpy, J/kg 

incremental compressor enthalpy, J/kg 

constant of integration 

kinetic energy, J/kg 

blade length, m 

Mach number 

specific speed 

total pressure ratio 

pressure, N/m 2 

volume flow rate, m 2 /sec 

dimensionless enthalpy (eq. (B75)) 

gas constant, J/(kg)(K) 

radius ratio 

Reynolds number 

radius, m 

temperature, K 



u 

blade speed, m/sec 

V 

absolute gas velocity, m/sec 

W 

relative gas velocity, m/sec 

w 

mass flow rate, kg/sec 

Z 

number of blades 

a 

absolute flow angle, deg from meridional 

& 

relative flow angle, deg from meridional 

% 

blade angle, deg from meridional 

y 

specific heat ratio 

6 

deviation angle between flow and blade, deg 

A 6* 

incremental boundary layer displacement thickness, r 

€ lim 

limiting impeller diameter ratio for slip calculations 

5 

vaneless diffuser loss coefficient 

V 

efficiency 

Atj 

decrement in efficiency 

X 

inducer hub-tip diameter ratio 

M 

2 

dynamic viscosity, (N)(sec)/m 

V 

2 

kinematic viscosity, m /sec 

P 

g 

gas density, kg/m 

fp 

flow coefficient 

'P 

head coefficient 

<j) 

angular velocity, sec 

Subscripts: 

AD 

adiabatic 

act 

actual 

aero 

aerodynamic 

av 

average 

BL 

blade loading 

cr 

critical state 

DF 

disk friction 
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est estimated 

H hub 

hyd hydraulic 

id ideal 

IGV inlet guide vane 
m meridional 

MF rms 

R rotor 

RC recirculation 
SF skin friction 

SL slip 

st stagger 

T tip 

th theoretical 

u tangential 

VD vaned diffuser 

VLD vaneless diffuser 

0 station just upstream of inlet guide vanes 

1 impeller inlet 

2 impeller exit 

3 vaned diffuser inlet 

4 vaned diffuser exit 
Superscripts: 

' absolute stagnation 

" relative stagnation 



APPENDIX B 


EQUATIONS 


The following are the equations listed in the order solved in the FORTRAN program: 

2ui 


a) = 


‘IT 


D 


IT 


(Bl) 


V cr = 


\M r 7 RT 0 

f y + 1 


(B2) 


V 1 rp 

V. = 1 V 

v IT v v cr 

v cr 


(B3) 


— sin 


u\ U 1T 


u/l V lT 


(B4) 


V mlT “ V 1T C0S ®1T 


(B5) 


W ulT = U lT 


1 - 



(B6) 


/3 1t = tan 


■l/ W ulT^ 

VmlTy 


(B7) 


W 1T ’ V^ulT + V mlT 


K = V mlT + 2V ulT = constant 


(B8) 

(B9) 


'W, 


W 2 = 


W 


tW 


IT 


it; 


(Bio) 


30 




(B11) 

(B12) 

(B13) 

(B14) 

(B15) 

(B16) 

(B17) 

(B18) 

(B19) 

(B20) 

(B2D 

(B22) 
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V s 


V mlMF = V K " 2V ulMF 


a 


IMF 


= tan 


-l/ V ulMF \ 
V mlMF/ 


V 2 2 

V ulMF + V mlMF 


W ulMF _ U 1MF ' V ulMF 


W 1MF = V V mlMF + W ulMF 


^ 1 MF = tan 


-1/ W ulMF ’ 


k V. 


mlMF/ 


T T' VlMF 

t imf- a o 2C 


^1T + ^1MF + ^1H 


5 lav 


-y/(y-l) 


n) 

\ p / IMF 


= 1 + 


IMF 


2C P T iMFy 


P 1MF ” P 1MF( 

\ p / IMF 


P 1MF " 


P 1MF 

RTlx 


p lMF " p 1MFL» 


1/y 


(B23) 

(B24) 

(B25) 

(B26) 

(B27) 

(B28) 

(B29) 

(B30) 

(B31) 

(B32) 

(B33) 

(B34) 


IMF 


(B35) 


Q = - d j t ( 1 + 


mlMF 


w - 


p 1MF Q 


. = _£o_ 

° RT^ 


Re = 


u 2 D 2 

MqPq 


rrM _ m . W 1MF 

1 =T imf + -^ 


2 2 

rr'* _ rrM , U 2 " U 1MF 

l 2 1 1 + 

2C„ 


W 2 

T 0 = T«’ 

z ^ 2C 


= 


<p = 


mlMF 


/w\ ,jj2 + pMF^ 
\ u 2/av 2 1 V D ' 


'W fl 


2 r 


\w 


IT/ 


2 / D 1T^ 

+ 

V Do I 


When inlet swirl is prescribed, equations (B45) to (B60) are solved instead of 
tions (B31) to (B44). 


a 


°st = 


IMF 


(B36) 

(B37) 

(B38) 

(B39) 

(B40) 

(B41) 

(B42) 

(B43) 

(B44) 

equa- 

(B45) 
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0.0076 


(B46) 


cos a 


1 + 


st ' 


^cos “ 0-025^ 


0.7 


KE = 


IMF 


“id" 


KE 
1 - e. 


£ 

\ p / IMF 


[i - _K^1 

L c p T oJ 


y/(r- 1) 


PlMT . L . 


Po 


C T« 

p U* 


J 1MF\ 


p’ -P’ VP ° 
P 1MF P 0 


A 

> p V IMF 


P 1MF 


P 1MF 

RT' 


~ 1 ( P \ 

P 1MF " p 1MF(J 

V p / IMF 


. /p\ 1/y 

p imf " p iMFrr 

V p / IMF 


w 


= PjQ 


( pV )lMF ” p lMF V m IMF 


(B47) 

(B48) 

(B49) 

(B50) 

(B51) 

(B52) 

(B53) 

(B54) 

(B55) 

(B56) 



(B57) 


J?L 


ip*v, 


cr/lMF 


y + l lv. 


2 

cr/O 


|l/(r-l) 


cr/O 


V n = 


cr/0 


cr 


M = 


9. 796X10" 7 
T 1MF + 110> 4 


Re = 


w 


/iD 


IT 


Ah 


0. 4 sin cy 


IGV 


IMF 


2 C0S ®1MF Re 


0.2 


( v o * v?mf) 


1+1 


IMF) 


Z = 6. 5 


D 


2 cosf^LlM 


1 -i 


''D 


IMF) 


Dr 


e lim " e 


-(8. 16 cos /3 2b /Z) 



(B58) 

(B59) 

(B60) 

(B61) 

(B62) 

(B63) 

(B64) 


e lim’ ^ en sR P velocity is expressed as 



V SL " “2 


^cos /3. 


1 - 


2b 


'IMF 

Do 


liml 


,0.7 


( 1 " e lim^ 


+ \lr 




IMF 

Do 


' e lim 


V ~ € lim ) 


. -l/ V SL C0S ^2b\ 
5 = sin 

l Wo / 


W 2id = W 2 


cos(|3 2b + 6) 


cos /3 


2b 


V u2 " u 2 " V SL " W 2id sin ^2b 


V m2 - W 2id cos ^2b 


= tan 


-1/ V u2 


m2/ 


W u2 = u 2 - V u2 


V 2 ’ V^m2 + V u2 


t 2 * T 2 + 


2C V 


ih aero = C p T i(^7 ‘ j 


Ah 


*th 


aero 


(B65) 

(B66) 

(B67) 

(B68) 

(B69) 

(B70) 

(B71) 

(B72) 

(B73) 

(B74) 

(B75) 
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_^2_ + 

0. 75 q th 

Wj T + W 2 

-/l - ° 1T ^ + 2 ° 1T 

Wit 

f V D 2 / D 2. 


Ah BL = 0.05D f 2 u| 


(B76) 


(B77) 


— 0. 02 DjUg ytan 

1 _ d imf 

L _ 0. 3048 

D 2 2 cos % 


(B78) 


(B79) 


Using the isentropic density rise as the first approximation, we solve equations 
(B80) to (B90) iteratively for the impeller exit density. 


p 2est " P 1MF( 


t_ \ V ( y - 1 ) 

— 

IMF/ 


Ah 


DF 


_»• 01356 p 2es t4 D l 


w Re 


0.2 


(B80) 


(B81) 


7rD 2 V m2 p 2est 


(B82) 


'IT 


D 


hyd _ 


D, 


D, 


2Z 


7T COS /3 


2b 


>2 1 - X 7r( 1 + X) 




+ X ‘ tan 2 /3 1T 


(B83) 


37 


JL_ 


ihsr=5 .6 C[ Dz /ty u| 

SF D hyd \2/ av 

(B84) 

D 2 


Ah act = Ah aero + Ah DF + Ah RC 

(B85) 

h id = Ah aero " Ah BL " Ah SF ' Ah IGV 

(B86) 

h id 

” B "ih 

aero 

(B87) 


p 2 




Ah 


^/(y- 1 ) 


aero 


C p T 0 


+ 1 ’ 


P 1MF 



M 


2 



r 2 



9. 7965X10' 7 T^ 5 
T 2 + 110.4 


(B88) 


(B89) 


(B90) 


(B91) 


(B92) 


(B93) 


- JL 
v - JZ - 

P 2 


(B94) 



(B95) 


AM 2 = 


A tan 


B 0 = L0 


F(l) = 


M 2 R 


1 + 


y - l 


m: 


o\y/(y _1 ) 

V 


C = 


C f r 2 

b«. 


A5* = 


/o. 37 r 2 AR \ 
y 8 cos a / 

* 2 r 2 AR \°- 2 


y COS O' 


b = b q 


2 A5 
bo 


AB = B Q - B 


(B96) 

(B97) 

(B98) 

(B99) 

(B100) 


- 2 ( 1+ ^M 2 ) 


(M 2 - sec 2 «) L 


(yM 2 - tan 2 o) S_ 


Bq cos a Bq AR 


?c 2 g j 

R J 


1 AB sec a M 2 AR 


(B 101) 


v = 


sec 2 ^ 


(M 2 - sec 2 cv) 


l + 2Lli M ^ S J-AB.^ltanaiR 

Bq cos a Bq AR R J 


(B102) 


El = 


P2 


1 



b 2 M 2 cos a 2 \ 2 



y/(y-l) 


(B103) 
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Equations (B98) to (B103) are solved by the numerical technique described in refer- 
ence 6 with increments in radius ratio of 0. 02 until a Mach number of 0. 8 or less is at- 
tained. 


p 3 


p 3 




(B104) 


Ah VLD = C p T 2 


ft> 3 \ (y ‘ 1)/y /P 3 \ (r " 1)/r ' 

W W 


(B 105) 


B t = 1 - B 


T I 

P 4 - P 3 


(B106) 

(B107) 


AH = 


M, 


r l + 


y - 1 


0.04' 


k (y+l)/2(y-l) 


P 4 


(B108) 


02 V^“i / 


ARNUM = p 4 AR 

(B109) 

P 4 = C p* (p 3 ’ p 3 ) + P 3 

(B110) 

/ v - 1 \rAr- 1 ) 

= P 4 ( 1+y 2 °- 04 ) 

(Bill) 

AR = ARNUM 

(B112) 


Equations (B108) to (B112) are iterated until successive approximations of area 
ratio agree within 0. 001. 
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Ah VD " C p T 2 


^\ (y-1)/y /PA (y ' 1)/y ' 


77 AD 


_ h id " Ah VLD " Ah VD 


Ah 


act 


PR = 


P 4 

PO 


H id = 77 AD Ah aero 


4 


N, 




Hr 


3/4 


id 


A 77 


Ah. 


BL 


BL 


Ah. 


act 


A?7 SF " 


Ah 


SF 


Ah 


act 


_ Ah DF 

At7 DF " 77 

Ah act 


a irc- RC 


Ah, 


‘act 


i„ - “VLD 
”VLD — 


‘act 


(B113) 

(B114) 

(B115) 
(B 1 16) 

(B118) 

(B119) 

(B120) 

(B121) 

(B122) 

(B123) 
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(B124) 


A? ?VD 

a? ?igv 


Ah 


VD 


Ah, 


'act 


Ah 


IGV 


Ah 


act 


(B125) 
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